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Abstract 

A cosmological model with a new variant of Chaplygin gas obeying an 
equation of state(EoS), P = Ap — where B = BqcP is investigated in 
the context of its thermodynamical behaviour. Here Bq and n are constants 
and a is the scale factor. We show that the equation of state of this ‘Vari¬ 
able Modified Chaplygin gas’ (VMCG) can describe the current accelerated 
expansion of the universe. Following standard thermodynamical criteria we 
mainly discuss the classical thermodynamical stability of the model and find 
that the new parameter, n introduced in VMCG plays a crucial role in deter¬ 
mining the stability considerations and should always be negative. We further 
observe that although the earlier model of Lu explains many of the current 
observational findings of different probes it fails the desirable tests of ther¬ 
modynamical stability. We also note that for n < 0 our model points to a 
phantom type of expansion which, however, is found to be compatible with 
current SNe la observations and CMB anisotropy measurements. Further the 
third law of thermodynamics is obeyed in our case. Our model is very general 
in the sense that many of earlier works in this field may be obtained as a 
special case of our solution. An interesting point to note is that the model 
also apparently suggests a smooth transition from the big bang to the big rip 
in its whole evaluation process. 
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1. Introduction 

Following the high redshift supernovae data in the last decade era we know that 
when interpreted within the framework of the standard FRW type of universe (ho¬ 
mogeneous and isotropic) we are left with the only alternative that the universe 
is now going through an accelerated expansion with baryonic matter contributing 
only 5% of the total budget. Later data from CMBR studies [2] further corrobo¬ 
rate this conclusion which has led a vast chunk of cosmology community ( j3] and 
references therein) to embark on a quest to explain the cause of the acceleration. 
In fact the studies on accelerated expansion and its possible interpretations from 
different angles have been reigning the research paradigm for the last few decades. 
The teething problem now confronting researchers in this field is the identification 
of the mechanism that triggered the late inflation. But as they are already discussed 
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extensively in the literature (we are sparing the readers here to repeat once again 
those arguments) all the alternatives face serious theoretical problems in some form 
or other. 

In the literature a good number of dark energy models are proposed but little is 
precisely known about its origin. Nowadays, the dark energy problem remains one 
of the major unsolved problems of theoretical physics [3]. On the way of searching 
for possible solutions of this problem various models are explored during last few 
decades referring to e.g. new exotic forms of matter ( e.g. quintessence) j5], phan¬ 
tom [6], holographic models [7], string theory landscape [8], Born-Infeld quantum 
condensate [§], modified gravity approaches ra. inhomogeneous spacetime HU, 
higher dimensional space time [12] etc. 

While the above mentioned alternatives to address the observed acceleration of 
the current phase have both positive and negative aspects a number of papers have 
come up taking into account the Chaplygin gas 1T51 flo] as a new form of matter 
held to simulate unified dark energy and dark matter model. It is presumed that 
for gravitational attraction, the dark matter component is responsible for galaxy 
structure formation while dark energy provides necessary repulsive force for current 
accelerated expansion. 

Motivated by the desire to explain away the observational fallouts better and 
better the form of the equation of state (in short, EoS) of matter is later generalised 
in stages first through the addition of an arbitrary constant with an exponent over 
the mass density, generally referred to as generalised Chaplygin gas (GCG) [3j. 
Barring the serious disqualification e.g., it violates the time honoured principle of 
energy conditions, it is fairly successful to interpret the observational results coming 
out of gravitational lensing or recent CMBR and SNe data in varied cosmic probes 
via the fine tuning of the value of the newly introduced arbitrary constant. The form 
of EoS is again modified through the addition of an ordinary mater held, which is 
termed in the literature as modified Chaplygin gas (MCG) [161 - 118] . claiming an even 
better match with observational results. 

It may be appropriate at this stage to call attention to a recent work by Guo and 
Zhang PM] where a very generalised form of the Chaplygin gas relation is invoked, 
assuming the constant B to depend on the scale factor a. Taking , e.g.,B = Boa~ n , 
it is shown that for a very large value of the scale factor the model interpolates 
between a dust-dominated phase and a quiessence phase (i.e., dark energy with a 
constant equation of state ) [I21H1] given by W = — 1 + |. 

As mentioned earlier, although different variants of dark energy models as well 
as Chaplygin type of gas models are brought in the cosmological arena as also 
its associated success to explain the observations coming out of different cosmic 
probes it has not escaped our notice that scant attention has been paid so far 
to address the important issue if all the so called perfect gas models are atleast 
thermodynamically stable. Otherwise they would lose the claim to be treated as 
a physically realistic system. Following this there has been of late a resurgence 
of interests among workers to address their queries to this aspect of the problem. 
Recently Santos et al have studied the thermodynamical stability in generalised [21] 
and modified Chaplygin [22] gas model on the basis of standard prescription [23j 
where both (i) (§p) 5 < 0, (§p) r < 0 and (ii) Cy > 0 are satisfied simultaneously. 

In the present work we attempt to generalise their results to Gou-Zhang formal- 



ism as also to clarify the issues concerning thermodynamics of the Chaplygin gas. 
We also investigate the nature of physical parameters such as pressure (P), effective 
equation of state (W), deceleration parameter (q) etc. on the basis of thermody¬ 
namics. Later many thermal quantities are derived as functions of temperature or 
volume. In this case, we also show as consistency check that the third law of ther¬ 
modynamics is satisfied by the new form of the Chaplygin gas. For the generalised 
Chaplygin gas we expect to have almost similar behaviour as the Chaplygin gas 
equations did show. Further, we see that the Chaplygin gas which shows a unified 
picture of dark matter and energy cools down through the adiabatic expansion of 
the universe without any critical point. Returning to the stability criterion of the 
Chaplygin gas we find a striking difference from the analysis of Santos et al. In 
our case of variable modified Chaplygin gas (VMCG) [23] we find that the stability 
depends critically on the new parameter n, as introduced by Guo et al. We also 
notice that stability decreases with the increase of magnitude of n, which apparently 
disfavors the formalism of Guo-Zhang [IP]. In the section on acoustic velocity we 
interestingly note that unlike the previous Chaplygin gas models where the squared 
sound velocity is always positive definite, the velocity here depends on the value of 
the newly introduced parameter n. But as noted earlier the stability criteria dictates 
that n should be negative which, however, makes the squared velocity also negative 
at the late stage of evolution. This means that the perfect fluid model for Guo’s 
variant of the chaplygin gas is classically unstable and is in line with the results 
obtained earlier by Myung [25] while dealing with the holographic interpretation 
for Chaplygin gas and tachyon. The paper ends with a short discussion. 


2. Formalism 

The line element corresponding to spatially flat FRW spacetime is given by 

ds 2 = dt 2 — a 2 (t ) ( dr 2 + r 2 d0 2 + r 2 sin 2 dd(j) 2 ) , (1) 

where a(t) is the scale factor. In this work we consider the following equation of 
state 

p = a p -A p) 

P a 

A and B are positive constants. As discussed in the introduction we have taken 
B = B 0 V~ a where n is an arbitrary constant and B 0 an absolute constant. Here P 
corresponds to the pressure and p the energy density of that fluid such that 

(3) 

where U and V are the internal energy and volume filled by the fluid respectively. 

We try to find out the energy density U and pressure P of Variable Modified 
Chaplygin gas as a function of its entropy S and volume V. From general thermo¬ 
dynamics, one has the following relationship 

= -P. 
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With the help of the above equations we get 



which, on integration yields 
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The parameter c is the integration constant, which may be a universal constant or 
a function of entropy S only. Now we rewrite the above equation in the following 
form 
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which has a dimension of volume. Now the energy density p of the VMCG comes 
out to be 
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From what has been discussed above we like to obtain the expression of relevant 
physical quantities and investigate their behaviour. 


(a) Pressure : 

Using equations (J2]) and (j9bj) the pressure P of the VMCG may also be deter¬ 
mined as a function of entropy S and volume V in the following form 
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The equation (HU gives a very general expression of pressure. Under special condi¬ 
tions 





































i) For n = 0 & A = 0, the equation (TTOT) reduces to 
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which is, as expected, similar to generalised Chaplygin gas (GCG) model [2Tj . 


ii) For n = 0&Aj^0,we get the modified Chaplygin gas (MCG) model with 
pressure given by m 
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iii) For + = 0, but n / 0, the model represents the Variable Chaplygin gas 
(VCGC) [26] and equation (TTOT) becomes 
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which is 


Again if we put a = 1 in equation (T13T) . we get P = — — r — ; —-tvt 

_ 2 { 1 + (f) }. 

identical with our previous work [27] when we consider the Variable Chaplygin gas 
(VCG) model. 

The pressure may have positive or negative values, depending on the magnitude 
of both A & n and also on volume V (fig - 1). For P — 0, let V is denoted by V c , 
which is given by 


V r = e 


3,4(1 +a) 
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(14) 


which restricts n as n < 3(1 + a). Initially, he., V < V c , P is positive, which indi¬ 
cates a radiation dominated universe. For V = V Cl P = 0 and V > V c , P is negative 
pointing to a state of accelerating universe. This is an interesting result showing that 
V c introduces a new scale in the analysis, beyond which a dust dominated universe 
enters the acceleration era. With expansion an initially decelerating universe tends 
to reverse its motion and prepares to accelerate when its volume crosses a critical 
value designated by V c . It is also to be noted that for physically realistic values of 
constants both V c and e are of the same order of magnitude i.e. e also signifies a 
volume scale beyond which accelerating era commences. SoF>e represents a very 
large volume and V <C e the reverse. In what follows we will see that this statement 
has significant cosmological implications. 






















(a) The graphs clearly show that pres¬ 
sure P can be both positive and nega¬ 
tive depending on the value of A. For 
A = 0, P is always negative, i.e. Chap¬ 
lygin type gas with n which is also in 
accord with the equation (|2j). 



(b) It show that pressure P is positive 
for small value of V and negative for 
large V when A > 0. As the value of 
n tends towards negative, P becomes 
more negative. 


(c) For A = 0, n = 0, P is always nega¬ 
tive, chaplygin type case. 


Figure 1: The variation of P and V for different values of A and n. Here we have 
taken B 0 — 1, a — 1 & c — 1 for constant S. 


(b) Caloric EoS: 


Now using the expression (|9bj) and (1T0)) we get the caloric equation of state 
parameter 


W = 


P 

P 



(15) 


As the last EoS is very involved in nature it is very difficult to extract much 
physics out of it. So we look forward to its extremal cases as: 

1. For small volume, V -C e, we get from equation fjl5|l that 


P ~ Ap. (16) 

This is a barotropic equation of state. In this case no influence of n on small 
volume. 

2. For large volume, V e, the equation (TT5|) reduces to 
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Equation (TT4l) shows that n < 3(1 + a). So there are three possibilities for >V 
depending on the signature of n as (i) n > 0, W > —1, here the caloric equation 
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Figure 2: The variation of TV and V for different values of n. The values of constants 

are taken as A = 1, Bq = l, a = l, c = 1. 


of state results in a quiescence type and big rip is avoided in this case, (ii) n = 0, 
W = —1, i.e., we get ACDM. (iii) n < 0, W < — 1 represents phantom like universe. 
This is compatible with recent observational results [28]. Influence of n is prominent 
in this case. 

Initially, he., when volume is sufficiently small, W > 0. As V increases to V c , 

W tends to zero. When V — V c — e , W = 0. Again V increases and 

W becomes negative. We have seen from hg-2 that n = 0, W = —1, i.e., ACDM 
model which is currently fashionable. But for n > 0 we get 0 > W > —1 for large 
V. Similar result was shown earlier in higher dimensional case also [l!2]. We do 


not discuss much about n > 0 because later we will show from thermodynamical 
stability considerations that the value of n < 0. 


(c) Deceleration parameter: 


Now using equation (fT51) we calculate the deceleration parameter of the VMCG 
fluid as 
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For mathematical simplicity here also we discuss the extreme cases. 
1. For small volume, V -C e, it gives 


1 3 . 
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1. e., q is positive, universe decelerates for small V. 

2. For large volume, V e, the equation (fT8]l reduces to 
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Initially, i.e., when volume is very small there is no effect of n on q, here q is positive, 
universe decelerates. But for large volume, q is negative and this depends on the 
value of n also. For flip in velocity to occur the flip volume (V/) becomes 






















Figure 3: The variation of q and V for different values of n. We have considered here 
A = ^, Bo = 1, a = 1 & c = 1. 
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Therefore n < 2(1 + a), which interestingly does not violate our previous restriction 
on n. A little analysis of the above equation shows that for Vf to have real value 
n < 2(1 + a), otherwise there will be no flip. This also follows from the fig - 3 where 
only n < 4 allows flip (for a = 1). Alternatively the inequality n < 2(1 + a) may 
be treated as our acceleration condition. Thus V < Vf, we get deceleration; and 
V > Vf , acceleration occur. 

It has not also escaped our notice that we are here getting two scales ( V c and Vf) 
for volume where apparently acceleration flip occurs. Again we know that in a FRW 
cosmology for flip to occur pressure should not only be negative but its magnitude 
should be less than | (i.e., p + 3 P < 0 ) as well. So in our case one should get 
V c < Vf , which also follows from their expressions (ITT)) and (I2TT) . In this connection 
it may also be mentioned that with V — Vf, p + 3 P = 0 as expected. 


(d) Velocity of Sound: 


Let us consider v s be the velocity of sound, then using equation (ITOl) we can write 
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Since sound speed should be 0 < (jffj < 1? now our analysis shows , 

1. For small volume, i.e., at early universe, the equation (1221) leads to 0 < A < 1, 
and this limit includes A — |, the radiation dominated universe. 

2. For large volume, the equation (1221) reduces to 
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The equation (1231) does not depend on A, it depends on n and a only. In what 
follows we shall see that from the thermodynsmical stability condition the value of 
n < 0, leading to a phantom universe [28J. Moreover the equation (T23l) gives an 
imaginary speed of sound for a > 0, leading to a perturbative cosmology. One need 
not be too sceptic about it because it favours structure formation }2TTj . 

It may not be out of place to draw some correspondence to a recent work by 
Y.S. Myung [25] where a comparison is made between holographic dark energy, 
Chaplygin gas, and tachyon model with constant potential. We know that their 
squared speeds are crucially important to determine the stability of perturbations. 
They found that the squared speed for holographic dark energy is always negative 
when imposing the future event horizon as the IR cutoff, while those for original 
Chaplygin gas and tachyon are always non-negative. This is in sharp contrast to our 
variant (VMCG) of the Chaplygin gas model where we observe that depending on 
the signature of n the squared velocity may be both positive or negative. However as 
discussed earlier a non negative value of n is clearly ruled out from thermodynamic 
stability considerations. This points to the fact that the perfect fluid model for 
VMCG dark energy model is classically unstable like the holographic interpretation 
for Chaplygin gas and tachyon and hence problematic in the long run. 

(e) Thermodynamical Stability: 

To verify the thermodynamic stability conditions of a fluid along its evolution, 
it is necessary (a) to determine if the pressure reduces both for an adiabatic and 
isothermal expansion [23] (§f) s < 0 & (§f) T < 0 and (b) also to examine if the 
thermal capacity at constant volume, cy > 0. 

Using equations (J2J) and ([XU we get 
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Now we have to examine the negativity of (fp) 5 - As this expression (j24l) is so 
involved we can not make any conclusion considering this equation as a whole. Our 
analyses naturally split into two parts : 


1. Firstly, we have considered small volume, V -C e, where equation (1241) gives 
(ff) ~ *-(l + A)y. We get from previous analysis that at the early stage 
of evolution P Ap , therefore, in this case (fp) 5 ~ — A(1 + A)^ which is 
independent of n but very much depends on A and at this stage of evolution 



2. For large volume, V e, the equation (1241) reduces to (|^) s , ~ ~ 3 FRT 0 )' 
Since P is negative at the late stage of evolution, so n must be negative to 
make ( j^) s < 0. We have seen that the dependence of n is prominent at the 
later case. From fig - 4 we get the similar type of conclusion. 


Now we discuss some special cases to constrain the parameters used here. 
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Figure 4: The variation of (§fJ) s and V for different values of n. The nature of graphs 
shows that for n < 0, (f^) 5 < 0 throughout the evolution unlike n > 0 where (§^) 5 < 0 
only at the early stage. We have taken A = 1 ,Bq = 1, a = 1, c = 1. 


(i) The simultaneous conditions n = 0, a = 0 and A — 0 must be discarded 
because it will place a severe restriction on the stability of this fluid, in such a case 
iw)s = ^ an d the pressure will remain the same through any adiabatic change 
of volume. However B 0 here behaves like a Cosmological Constant. We get the 
de-Sitter type of metric. 

(ii) Again for a = 0, A = 0 and n 0, (fp) s = § B 0 V~^ 1+ %\ i.e., for n < 0, 
(w) S < ^ which is i n disagreement with the previous work of Santos et al [21 j 
where simultaneously a = 0, A = 0 is not possible. Relevant to mention that, n < 0 
implies that the pressure goes more and more negative with volume. This agrees 
with the observational results [19], 

(iii) When A = 0, n = 0 and a 0, the equation (124|) reduces to Santos’s work 
[21j. In this case (§f) s will be 
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which gives (§p) 5 < 0 for a > 0. ft also agrees with the work of Santos et al [21J 
in this held. 

(iv) when a = 0, n = 0 and Aff 0, (f^) s = — (y) , i.e., (§f) s < 0 which 

implies simultaneously A > 0 and Bq > 0. 

(v) For A 0 but a 0 & n = 0 the equation (1241) reduces to 
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This equation (126|) is identical with the another work of Santos et al 

(vi) Again when A = 0 but the equation (1241) gives 
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This is the case of Variable Generalised Chaplygin gas (VGCG) model [[26] . 
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Figure 5: The variation of s and V is shown for A = 0, a = 1.53, n = 0.75 and 

also Bq = 1, c = 1. 

It is clear from equation f|27|) that for a > 0, n should be negative for (|^) < 0 
throughout the evolution. It may not be out of place to call attention to an earlier 
work of Lu [26] where he studied the Union SNe la data and Sloan Digital Sky Sur- 
vey(SDSS) baryon acoustic peak to constrain the Variable Generalised Chaplygin 
Gas (VGCG) model where he obtained the best fit values of n — 0.75 and a = 1.53. 
Our analysis shows that Lu’s conclusion is untenable if consideration of thermody¬ 
namical stability is taken into account. Similar conclusion may be drawn from the 
analysis of our fig-5. In this connection it may be pointed out that we have studied 
pza for the case of A = 0 but a = 1 & n ^ 0, where we have seen that (§f) s < 0 
for n < 0 throughout the evolution. 

But our analysis is based on all these constants: A, Bo, a and n. As it is not 
possible to constrain all the parameters simultaneously from equation (l24jh we have 
tried it step by step. Now we want to concentrate on the signature of n rather 
than the other parameters because we are dealing with VMCG model where n has 
a special role throughout the evolution, particularly at late time. Now n = 0 was 
discussed earlier by several authors [22]. It is seen from the equation (124|) that for 
positive values of A, B 0 and a, it will be (f^) 5 < 0 when n < 0 throughout the 
evolution. Fig-4 gives similar type of conclusion. Again it is seen from the fig-4 that 
for the positive value of n at small volume (§p) 5 < 0 but as volume increases, i.e., 
at the late universe where the influence of n is significant, for n > 0, (§f^) 5 > 0; 
so the stability is questionable throughout the evolution for n > 0. Therefore, in 
the context of thermodynamical stability, we have to conclude that the value of n 
should be negative for the VMCG model. 

Now we have to examine if (§f?) r < 0 as well. We will show in the next section 
that for n < 0 this condition may also be satisfied. 

One should also verify the positivity of thermal capacity at constant volume c„ 
where c v = T (ff) r = (f^) r = V (§r)y Now we determine the temperature T of 
the Variable modified Chaplygin gas as a function of it volume V and its entropy S. 
The temperature T of this fluid is determined from the relation T = (§^) r - Using 
the above relation of the temperature and with the help of equation (]6]) we get the 
expression of T as 
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If c is also assumed to be a universal constant, then = 0 and the fluid, in 
such a condition, remains at zero temperature for any value of its volume and pres¬ 
sure. Therefore, to discuss extensively the thermodynamic stability of the variable 
modified Chaplygin gas whose temperature varies during its expansion, it is neces¬ 
sary to assume that the derivative of equation (I25|) is not zero implying (J|) ^ 0. 
We have no a priori knowledge of the functional dependence of c. From physical 
considerations, however, we know that this function must be such as to give posi¬ 
tive temperature and cooling along an adiabatic expansion and so we choose that 

(I) > 0- 

Now from dimensional analysis, we observe from equation (J5J) that 

M = {yjlU} 1+ "- < 29 > 

Since [U] = [ T} [S'], we can write 


[c] = [T] 1+a [S] 1+a [y] A(1+a) . (30) 

It is difficult to get an analytic solution of c from equation (13Uj) . so as a trial case, 
we take an empirical expression of c and then to check if the resulting expressions 
satisfy standard relations of thermodynamics. But as c is a function of entropy only, 
the expression of c will be 


c= (rv A ) 1+a S 1+a , (31) 

where r and v are constants having the dimensions of time and volume respectively. 
Now 
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Using equation (j28l) and (132]) . we get the expression of temperature 
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and from equation (I33a[) . the entropy is 
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It follows from equation (jM[) that for positive and finite entropy one should have 
0 < TV a < tv a , but individually 0 < T < r and v < V < oo, i.e., r represents the 
maximum temperature whereas v represents the minimum volume. It is shown from 
the equation (133cl) that as volume of the VMCG increases, temperature decreases. 
It is also proved from the equation (I33cj) that when T —> 0, V —> oo and when 
T —y t, V —> v. Thus we can apparently avoid the initial singularity. 

Evidently at T = 0, S = 0 which implies that the third law of thermodynamics 
is satisfied in this case. 

Now using equation (1341) we get the expression of thermal heat capacity as 



(35) 


Since 0 < TV A < tv a and a > 0, cy > 0 is always satisfied irrespective of the 


value of n. This ensures the positivity of a. It is interesting to note that when 
the temperature goes to zero cy goes to zero as expected from the third law of 
t her mo dynamics. 

If we put A = 0 and a = 1, i.e. Variable Chaplygin gas model, we get the 
identical expression of cy of our previous work [27]. Again for A = 0 and n — 0, 
the equation (]35|) reduces to the work of Santos et al [2TJ . 

To end the section a final remark may be in order. While positivity of specific 
heat is strongly desirable vis-a-vis when dealing with special relativity, in a recent 
communication Luongo et al [30j argued that in a FRW type of model like the 
one we are discussing a negative specific heat at constant volume and a vanishingly 
small specific heat at constant pressure (cp) are compatible with observational data. 
In fact they have derived the most general cosmological model which is agreeable 
with the cy < 0 and cp ~ 0 values obtained for the specific heats of the universe 
and showed, in addition, that it also overcomes the fine-tuning and the coincidence 
problems of the ACDM model. 

(f) Thermal EoS: 

Since P = P(V,T ), using (£>]), (13TT) and fISTjl we get the internal energy as a 
function of both V and T as follows: 



( 36 ) 


Now using (j2J), © and (155]) the Pressure is 






















p = p 


(37) 
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We have seen from equation ()37l) that for A = 0 and a = 1, the solution reduces 
to our previous work m It is to be seen that for A = 0 and n = 0 the above 
solution goes to an earlier work of Santos et al m 
Now the thermal EoS parameter is given by 



\ Qt 




(38) 


Thus thermal EoS parameter is, in general, a function of both temperature and 
volume. If we consider early stage of the universe when at very high temperature 
and small volume, i. e., at T —> r and also V —> v, we get from equation (j38j) that 
u) ~ A, he., P ~ Ap. This is same as equation 



0 1 2 3 4 5 

V -> 

Figure 6: The variations of (§y) T vs V are shown. The graphs clearly show that 

(§F )t < ^ throughout the evolution only for negative values of n. Here we have taken 
A = 0.1, a = 0.1, Bq = 1, r = 1 and v = 1. 

Secondly we consider for large volume, he., very low temperature , he., T —>■ 0, we 
get from equation (138|) that uj ~ —1 + , which is identical with equation ()2T]) as 

is customary with the existing literature in this held [snug. Thus thermodynamical 
state represented by equations (TT6|) and (fl7j) are essentially same at both early and 
the late stage of the universe. 

Now from equation (j37j) we have to examine if (|^) T < 0, and a lengthy but 
straight forward calculation shows that only for negative value of n, this condition 
satisfies. It is very difficult to find out (§p) r from the equation (1371) in a compact 
form and to get requisite inferences from it. So we have taken recourse to graphical 
method instead. From fig - 6, we clearly infer that (fp) T < 0 for n < 0 throughout 
the evolution depending on the values of constants. 

Now to discuss the thermodynamic stability of the VMCG one should have 
(§t)t < 0 anc ^ (w) s < *' e ’’ both isothermal and adiabatic situations need 

to be addressed. In the process we have found that (§p) 5 and (|^) T are negative 
for n < 0. Relevant to mention that in a previous work by Santos et al [22] it 
















was assumed (§f) T = 0 for simplicity, but in our case we have not made such a 
simplistic approach. In fact we have used a different type of approach to find out 
the expression of c = c(S). As we have not made this assumption our analysis is 
more general in nature and our formulations would not reduce to there case when 
n = 0. 

As pointed out in the previous section we here take up a standard relation of 
thermodynamics (e.g. the first internal energy equation) [31] as a trial case to see 
if the equation (13T[) is satisfied. 


fdU 

[dv 


T 


T 



(39) 


Using equations (1361) and (1371) we verify the relation (1391) which proves the cor¬ 
rectness of our approach. 


(f) Pressure-Volume relation: 


It is very difficult to arrive at a (P ~ V) relation in a general way. So one has 
to take recourse to extremal conditions. 

1. For small volume, f «e, which gives P ~ Ap. In this case energy density 
and also pressure are very high. Using (I33cl) the temperature will be 


T « 


tv a 

W' 


(40) 


At the early stage of the universe V —> v (minimum volume) implies that T —>■ r 
(maximum temperature) So the temperature is high enough at this stage. Using 
equations (I33cl) and (140]) we get 


therefore, 


P 


S 


tv a 

V 1+A 
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UV A = pV A+1 = Stv a . 


(41) 

(42) 


We know at this stage that P Ap = Ay, i.e., PV AU, which gives using 
equation (1421) PV 1+A = Srv A . Since the entropy remains constant in an adiabatic 
process, this relation leads to PV 1+A = Constant. So it is observed that for small 
volume, i.e., at high temperature the VMCG behaves as a fluid of y(= ^) = 1 + A. 
We can also rewrite the EoS as P — (7 — 1) p. Since early universe is radiation 
dominated i.e., A = ~, the value of 7 = l + A = |, the pressure is related to 

4 

the volume as PV 3 = constant. Thus the VMCG behaves like a photon gas. The 
equation of adiabatic photon gas coincides with extreme relativistic electron gas. 

2. For large volume, V » e, 

Due to low density at this stage, entropy density is sufficiently small at the low 
temperature. We know from equation (fT71) that 
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3(1 + a) 

Again, from equation (jflbjh in this case we get 
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Now using equations (j43l) and (J44j), we get 
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l + CK 


(44) 


py 3(i+ a ) — ) _i _)_ 
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B 0 (l + a) 1 1+ 


3(1 + a) J [ N J ' (45) 

Equations ()43lh (|44li and (|45li can be obtained from the equations (|38jl . ()36li and 
()37li respectively at T —>■ 0. For an adiabatic system, the above equation looks like 
PV 1 = k , where k is a constant. Since we know that VMCG is thermodynamically 
stable for n < 0 which leads to k < 0. In this case, P is also negative. At the 
late stage of evolution, i.e., at low temperature, the VMCG behaves like a fluid of 
7 — 3 ( (+o) • Interestingly, it is to be noted that the value of this 7 depends on both 
n and a at this stage of evolution. It further follows from equation (133(1) that T —> 0 


for y << 1. Since any system near T = 0 is in states very close to its ground state, 
quantum mechanics is essential to the understanding of its properties. Indeed, the 
degree of randomness at these low temperatures is so small that discrete quantum 
effects can be observed on a macroscopic scale. 


3. Discussion 

Following the discovery of late acceleration of the universe there is a proliferation 
of varied dark energy models as its possible rescuers. While many of them sig¬ 
nificantly explain the observational findings coming out of different cosmic probes 
serious considerations have not been directed so far to the query whether the models 
are thermodynamically viable, for example if they obey the time honoured stabil¬ 
ity criteria. We have here considered a very general type of exotic fluid, termed 
‘Variable Modified Chaplygin gas’ and studied its cosmological implications, mainly 
its thermodynamical stability. Regarding the cosmological dynamics we have come 
across two characteristic volumes of the fluid - V c and Vf representing critical volume 
and flip volume respectively. The former refers to the case where pressure changes 
its sign while the latter gives the volume when the acceleration flip occurs. From 
physical consideration one should get V c <Vf - which also matches with our analy¬ 
sis. As discussed at the end of section 2(c), for flip to occur in FRW cosmology only 
a negative pressure is not the necessary and sufficient condition. The magnitude of 
pressure should be also less than |, which also follows from our analysis. 

Although the exhaustive analysis of the latest cosmological observations provides 
a de fini te clue of the existence of dark energy in the universe but it is difficult to 
distinguish between the merits of various forms of dark energy at present. For the 
stability criteria we have followed the standard prescription : (fp) 5 < 0, (|^) T < 0 
and cy > 0. This, however, dictates that the new parameter, n introduced VMCG 
should be negative definite. This contrasts sharply with an earlier contention of 
Lu’s [26] where to explain the observational results they have to choose n as posi¬ 
tive definite. This is pathological because it makes the system thermodynamically 
unstable. 








Again this model shows that at early stage, the EoS becomes P = Ap, where 0 < 
A < 1. But at late stage, it reduces to equation (l43]h From the thermodynamical 
stability conditions, we find that n < 0, which favours a phantom like evolution and 
big rip is thus unavoidable. So far as phantom model is concerned, it is found to 
be compatible with SNe la observations and CMB anisotropy measurements [28] , 
The most important conclusion coming out of our analysis is that this model covers 
both big bang and big rip in the whole evolution process. 

It is to be noted that at T = 0 the entropy of VMCG vanishes as in conformity 
with the third law of thermodynamics. We have studied both the thermal and the 
caloric EoS which shows that both 0 < T < r and v < V < oo where r is a 
maximum temperature and v is a minimum volume attainable. Here r and v are 
canonical in the sense that as T —» r, V —> v at the early stage, i.e., r represents 
the maximum temperature that our VMCG model can sustain for the small finite 
volume, v. 

We have also discussed (P ~ V ) relation and at early stage (A = |) it is shown 
that for an isentropic system VMCG behaves like a photon gas, as it was at the 
time of radiation dominated era. It is also shown that for large volume, i.e., at 
low temperature, entropy is sufficiently low which agrees well with the currently 
available low energy density of the universe. 

Finally it may not be out of place to point out that as the field equations are 
very involved in nature we have to adopt an ansatz to determine the function of 
integration c = c(S). However to justify it we have checked an important relation 
(1391) (first internal energy equation) and have found it to be correct. So our ansatz 
is essentially viable. Moreover, this model is very general in the sense that many of 
earlier works in this field may be obtained as a special case. 

To end a final remark may be in order. We have here concentrated on the 
cosmological and thermodynamical behaviour of the VMCG model mainly on a 
theoretical premise. However, as a future exercise one should try to constrain the 
value of the parameters associated with both thermal and caloric EoS in the light 
of observational values. 
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